Nodal Quasiparticle Dispersion in Strongly Correlated d-wave Superconductors 
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We analyze the effects of a k-dependent self-energy on the photoemission momentum distribution 
curve (MDC) dispersion and lineshape. We illustrate this general analysis by a detailed examination 
of nodal quasiparticles in high Tc cuprates. Using variational results for the nodal quasiparticle Z, 
which varies rapidly with hole doping x, and v 1 ™ , which is independent of x, we show that the high 



energy dispersion v^ gb 



" /Z, so that it is much larger than the bare band structure dispersion 



and also exhibits strong doping dependence in good agreement with recent photoemission data. 
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Nodal quasiparticles (QP's) are the dominant low-lying 
excitations in the superconducting state of the high Tc 
superconductors, and have been the focus of intense ex- 
perimental and theoretical investigation I n particu- 
lar, angle resolved photoemission spectroscopy (ARPES), 
exploiting the simple Lorentzian lineshape of momentum 
distribution curves (MDC's) |2|, has provided a great 
deal of insight into the nodal QP weight and dispersion 
H H H Ull II H Two important questions are: (i) 
What constraints do these experiments place on a model 
for the superconducting state? (ii) To what extent are 
nodal QP properties determined purely by strong elec- 
tronic interactions? 

The purpose of this Letter is twofold. First, we exam- 
ine the effect of the k-dependence of the electronic self 
energy on the analysis of MDC's for gapless nodal QP's. 
Specifically, we determine how this affects the MDC dis- 
persion, after clarifying the constraints under which the 
MDC lineshape is Lorentzian. This is important, because 
the usual basis for understanding a Lorentzian lineshape 
is a k-independent self-energy which is certainly suf- 
ficient but by no means necessary as shown here. 

The second goal is to use the above formalism to gain 
insight into some unusual aspects of the nodal QP dis- 
persion in the high Tc superconductors. Some time ago 
we predicted, on the basis of a variational calculation[10j 
for a Hubbard model that the nodal QP spectral weight 
Z vanishes linearly with hole doping x even though 
their Fermi velocity v l p w remains doping independent as 
x — > 0. (We use the notation v l p w here, rather than 
the more conventional vf, to clearly distinguish it from 
v F lgh defined below). These predictions have now been 
verified in some detail by angle-resolved photoemission 
spectroscopy (ARPES) experiments 0,0, 0|- 

These same experiments also show an unusual doping 
dependence of the high energy dispersion, at bin ding en- 
ergies larger than the much-studied kink [3, la, Ll IllL LL2| 
in the nodal MDC dispersion. The "high energy" Fermi 
velocity w^ lgh , at energies of around 150 to 200 meV, is 
experimentally found 0, HQ to be quite strongly dop- 
ing dependent and increases markedly with underdoping, 
where as one might have naively expected the "bands" 



to become narrower with decreasing x. This behavior 
of i>^ Igh is also in marked contrast with both the Fermi 
velocity at Ep which, as stated above, is x- independent, 
and the "bare" band-structure Fermi velocity v F which 
has a rather weak x-dependence. 

In this Letter, we show that this surprising behavior of 



(x) follows directly from the same theory 10] which 



high 

predicted Z(x) and the x-independent v 1 ™ . As discussed 
below, the results for Z(x) and v l p w imply strong con- 
straints on both the u>- and k-dependence of the real part 
of the self energy E'(k, u>). Specifically, both dE'/du) and 
dT//dk must exhibit 1/x singularities. The strong k- 
dependence is somewhat unexpected, and we show here 
that it has an important effect on the "high energy" dis- 
persion: 
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dZ'/dk = v l ° w /Z 



(1) 



leading to a strong doping dependence in v F lg,h , including 
an 1/x divergence as the hole doping x — > 0. We make 
detailed quantitative comparison of our results with ex- 
isting ARPES data, and conclude with comments on the 
kink in the MDC dispersion. 

MDC analysis: We first develop the formalism for the 
analysis of photoemission MDC's, focusing on two points 
not emphasized in the literature: the conditions under 
which the MDC lineshape can be Lorentzian despite k- 
dependence of the self energy, and the role of dT,'/dk 
in the MDC dispersion. Recall that the ARPES in- 
tensity from a two-dimensional system is given by [l^j ] 
J(k, u>) = Io(k)/(w)Ak, w where Iq contains the dipole 
matrix element and kinematical factors, f{uf) is the Fermi 
function and A^ ^ is the one-particle spectral function. 
Here k is the 2D momentum and u> the energy measured 
from the chemical potential. The usual energy distribu- 
tion curve (EDC) method views ARPES data as a func- 
tion of energy at fixed k. The MDC method, on the other 
hand, analyzes /(k, uS) as a function of k, along a suitable 
chosen cut in k-space, for fixed u>. The MDC method 
is more powerful than the traditional EDC for study- 
ing gapless excitations, because of the simple Lorentzian 
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MDC lineshape and also the ease of unambiguously iden- 
tifying the extrinsic background. 

Here we are specifically interested in the MDC's for 



a dr r 



superconductor for k along the zone diagonal, 



so that the "off-diagonal" self energy (gap function) van- 
ishes. The formalism derived is obviously applicable near 
the Fermi surface of any normal metal also. The spectral 
function is then given by: 



|E"(k, W ) 



7T [w - e (k) + no - S'(k, uj)} 2 + [E"(k, uj)}< 



(2) 



where £o(k) and fio are the "bare" dispersion and chem- 
ical potential, and E = E' + iE" the self-energy. 

The "renormalized" fe F is defined by the condition 
e (k_F) — /-to + E'(kp , 0) = 0. In case interactions lead to 
no change in the fop, this places a constraint E'(k F , 0) = 
on the real part of the self energy. We now expand 
in small \k — kp\ <C kp in the vicinity of the node at 
fcp, keeping to arbitrary for the moment. Linearizing the 
dispersion, we thus get 



e (k) - + E'(k F , 0) ~ v° F (k - k F ) 



(3) 



where v F is called the bare Fermi velocity. 

We next assume E"(k, u>) ~ Y,"(k.F,td) which is equiv- 
alent to setting the first correction 



dT,"(k F ,cj)/dk = 0, 



(4) 



for all to of interest. This is a necessary condition for 
a Lorentzian lineshape in (k — kp), since without it the 
numerator in eq. J2J) would have a term linear in (k — kp) 
which would make the expression non-Lorentzian. The 
Kramcrs-Kronig transformation of the above condition 
implies that <9E' jdk can only be a constant independent 
of w, a condition we will need to use below. (In the 
experiments, an additional constraint for a Lorentzian 
MDC is that the matrix element prefactor Iq should be 
k-independent in the range of relevant k's.) 

With this one assumption, the spectral function Ak^ 
is easily shown to be a Lorentzian in (k — kp) peaked at 



k(u>) = k F + 



-[E'(k F ,w)-E'(k F ,0)] 



(v° F 



dE'/dk) 



(5) 



with a width Afc = \E"(k F ,cj)\/(v F + dT,'/dk). The 
dispersion of this Lorentzian MDC is then obtained from 



dk _ l-dV'(k F ,w)/dw = CM 

~cL, ~ v° F + dE'/dfc = vpTdWJdk' 



(G) 



where the denominator is w-independent and evaluated 
at k=k F , but the numerator at arbitrary lu and k=kp. 

In the limit oj — > 0, the MDC dispersion of eq. © yields 
the standard result. Defining the low-energy renormal- 
ized Fermi velocity v 1 ™ 
get 



via l/v l ° w = dk/duj{uj -> 0), we 
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(7) 



Here the QP weight Z is defined by 

Z = 1/C(w -> 0) = 1/ [1 - dE'/dw] 



(8) 



with 9E' /du) is evaluated at (k = k F , to = 0). 

On the other hand, one can quite generally show [y| 
that at intermediate to high energies, \dT,'/duj\ <C 1, 
so that w 1. Then the corresponding high energy 



l/i> F gh = dk/did is given by 



high 



dZ'/dk, 



(9) 



which together with eq. Q leads to the result eq. (|TJ 
stated in the Introduction. Here dYj'/dk is an id- 
independent constant as discussed in connection with 
the necessary conditions (Kramers-Kronig transform of 
eq. for Lorentzian lineshapes. 

The above results are completely general, and we next 
illustrate this formalism for the low and high energy dis- 
persion of nodal QPs in the cuprates. 

Low energy properties Z(x) and v l F w (x): In Ref. 01 
we have presented a T = variational theory, building on 
early ideas of Anderson for the strongly correlated 
d-wave SC and its low-lying excitations, examining how 
they evolve as a function of hole doping x from a Fermi 
liquid state at overdoping x £ 0.35 to a Mott insulator 
at half-filling (x = 0). Here we summarize specific re- 
sults from Ref. which relate to doping dependence 
of the nodal QP spectral weight Z(x) and Fermi velocity 
v l p w (x), and describe their implications for E'(k, oj). 

We begin by emphasizing that these results for the 
nodal QP were not obtained by working with a varia- 
tional excited state, e.g., a Gutzwiller projected Bogoli- 
ubov QP. Given the very broad linewidth of the spectral 
function, it is difficult to use simple variational excited 
states to obtain useful results on the coherent QP piece; 
see Ref. |16|. Instead we focused on the singularities as 
a function of k in the energy-integrated moments of the 
spectral function to extract information related to gap- 
less QPs. In particular, Z was extracted from the jump 
in cLuA^^ — n(k) and v l p w from the slope disconti- 
nuity of didioAvui at k = k F 

We found that [lfj Z(x) decreases with underdoping, 
with Z ~ x as x — ► 0, and remarkably that v l p w (x) 
is essentially doping independent and remains finite as 
x — > 0. Using eq. 10 one must then conclude from the 
calculated Z(x) that |<9E'/<9w| ~ 1/x as x — > 0. Further, 
using eq. (0 , the calculated v l p w (x) implies a compensat- 
ing divergence dT,'/dk ~ Ja/x, where the superexchange 
J and lattice constant a enter on dimensional grounds. 

In the x — * limit one thus obtains Z ~ x, with a 
renormalized low energy Fermi velocity v l p w (x — ► 0) = 
CJa, where C is a dimensionless constant of order unity. 
Remarkably these results are independent of the bare 
band-structure dispersion and thus "universal" . From 
our numerical results |lfj| we find C ~ 4.5. We emphasize 
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FIG. 1: (a) Doping dependence of the nodal QP weight Z in 
Bi2212 extracted from data in Ref.® using Z = 1/(1 + A). 
The doping x was calculated from their sample T c using the 
empirical formula 86(x - 0.16) 2 = (1 - T c /T™ ax ), with 
T ma X = Q1K for their samples ( b ) z ^ predicted from a 

variational approach shows that Z — » on underdoping, and 
is in good agreement with the experimental trend and values. 

that these non-trivial results for the doping variation of 
Z and v l ™ are non-perturbative, strong coupling results 
which come from a variational calculation which properly 
takes into account the strong local Coulomb repulsion. 

High Energy Dispersion vp lsh (x): Combining all the 
results described above, we see from v F lsh = v l p w /Z(x), 
that we expect vp lsh ^> v l ™ and also that it will show 
considerable doping dependence and increase with un- 
derdoping, given the predicted Z(x). Moreover, using 
eq. JjJJ and the form of <9£' /dk derived above, it immedi- 
ately follows that there is Ja/x divergence in v F lsh , which 
at small enough x is much larger than and independent 
of the "bare" band-structure v F . 

Comparison with ARPES Experiments: We are 

now in a position to compare the theoretical results de- 
rived above with existing ARPES data. In Fig. 1 we plot 
the nodal QP weight Z as a function of hole doping x. 
Fig. 1(a) shows Z extracted from the experimental data 
of Ref. [i| on Bi2212. These authors plot A = -dT,'/dco 
as a function of doping, from which we obtain Z using 
Z = 1/(1 + A). In Fig. 1(b) wc plot the nodal Z(x) ob- 
tained from our variational calculation an d find very 
good agreement between theory and experiment. Wc 
should note that to leading order in J/t [ITJ the calcu- 
lated Z(x) is independent of the choice of input parame- 
ters in our theory |2fJ. The nodal Z(x) has also been re- 
ported for LSCO but in arbitrary units, which makes 
a quantitative comparison difficult, however the trend is 
qualitatively similar to our theoretical prediction. 

We next show in Fig. 2 the low energy Fermi velocity 
v X p™ for nodal QPs in units of eV-A as a function of dop- 
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FIG. 2: (a) The Fermi velocity v F in Bi2212 (open squares 
joined by a line as guide to the eye, Ref.[(|) and LSCO (open 
triangles, Ref. 7°]) is nearly doping independent with a similar 
value in various cuprates (see RehpJ for additional data), (b) 
Theoretically predicted v 1 ™ using a variational approach is 
constant for x<0.2, independent of the bare band structure 
( "universal" ) , and crosses over to the bare Fermi velocity (v F , 
dashed line) at large x. The values and trend are in good 
agreement with experiment Q. 



ing. The experimental results shown in Fig. 2(a) are ob- 
tained from the MDC dispersion at Ep for LSCO (open 
triangles from Ref. 0) and for Bi2212 (open squares from 



Ref. In Fig. 2(b) we show the "bare" 



(dashed 



line), which has a weak doping dependence coming from 
the change of the chemical potential in the bare tight- 
binding model |20j . The renormalized low energy Fermi 
velocity v l F w is however found to be smaller and nearly 
doping independent, within the error bars of the theoret- 
ical calculation for the hole doping range x < 0.2, 
its scale being set by Ja as explained in detail above. 

Finally we come to the comparison of theory and ex- 
periment for the high energy dispersion v F lsh in Fig. 3. 
The data in Fig. 3(a) are obtained from the high energy 
MDC dispersion for LSCO (open squares from Ref. |7|) 
and for Bi2212 (open triangles from Ref. @). In Fig. 3(b) 
we plot the theoretical estimate of the high energy veloc- 
ity using the expression derived above iy lgh = jZ{x). 
We find good agreement with the data over the range of 
dopings from x > 0.05 both in terms of magnitude and 
overall trend. There is a 1/x divergence in the theoret- 
ical value at small x, for reasons explained above, but 
evidence for that in the very low doping data in LSCO is 
not clear. Quite independent of our variational calcula- 
tion, the general result vp lsh = vp™ /Z(x) derived above, 
together with the experimentally observed behavior of 



.low 



(constant 0) and of Z(x) (vanishing as 



8]) 



0. 



in LSCO, implies a w^ lgh that increases strongly as x 
We hope future experiments will clarify this situation. 



Implications for the kink: The success of the theory 
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FIG. 3: (a) Doping dependence of the high energy energy 
velocity t;^ lgh above the kink in the MDC dispersion from 
data on Bi2212 (open squares joined by a line as guide to 
the eye, Ref.®) and LSCO (open triangles, Ref.0). (b) A 
variational calculation |l9j| of w^ lsh = i>j? w )Z shows that it 
increases on underdoping in agreement with the experimental 
trend and values for x £ 0.05. 



described above leads to the obvious question: What is 
the origin and energy scale of the kink? By the "kink" 
we mean the sharp change in the observed MDC disper- 
sion separating the low and high energy regimes which 
occurs at an energy scale of about 70 meV ja, IE • Its 
origin is controversial at present, with electrons interact- 
ing with either a n op tical phonon [f| or with the neutron 
resonance mode Qjlllll being the two likely scenarios. 

We have presented here a framework for understanding 
the low and high energy limits MDC dispersion of eq. Jfjjl, 
in limits where the factor of C,(oj) goes to Z and unity 
respectively. This general framework cannot directly ad- 
dress precisely how and at what energy scale there is a 
transition from one limiting behavior to the other. The 
variational approach only permitted a calculation of the 
low energy behavior of the nodal QPs — its implica- 
tions for i^ Igh hinged upon the remarkable relationship 
between this high energy quantity and and the low en- 
ergy v l p w and Z, which could be calculated. However, 
we emphasize that the variational calculation, which fo- 
cuses on the strong Coulomb interactions, is able to give 
qualitative and semi-quantitative insights into the doping 
dependence of both the low and high energy dispersion 
of the nodal QPs. Thus it would be very surprising if the 
intermediate energy scale kink were not also dominated 
by strong electron-electron interaction effects. 

Conclusions: To summarize, we first obtained the con- 
ditions, eq. and its Kramers-Kronig transform, for a 
Lorentzian MDC lineshape, and then used this to derive 
the MDC dispersion eq. JHJ) and the linewidth, both of 
which involve the k-dependence of the self energy in an 
essential way. We next showed how the high energy dis- 



persion was related directly via eq. to the QP spec- 
tral weight Z and the low-energy Fermi velocity v l p w . 
We then illustrated this general formalism by using vari- 
ational results for the nodal QP in the cuprates, where 
we showed that the Z(x) was strongly doping dependent 



and vanishing as x — > 0, 
dependent, and Vp S = 



.low 



,low 



was essentially doping in- 
jZ[x). All of these results 



were shown to be in good agreement with ARPES data. 
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